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ABSTRACT 


This  report  presents  the  results  of  theoretical  and  experimental  investigations 
of  the  flow  in  the  wake  of  hypersonic  vehicles.  Nonsimilar  solutions  of  the 
equations  governing  the  "far  wake"  are  obtained  and  contrasted  with  similar 
solutions  and  with  solutions  involving  the  Karman-Pohlhausen  approach.  Nu¬ 
merical  examples  of  the  flow  field  in  the  wake  of  ballistic -range  and  wind- 
tvmnel  models  and  of  a  flight- size  vehicle  are  presented. 

A  series  of  ballistic -range  experiments  with  conical  models  of  10,  15,  and 
27-1 /2-degree  half-angles  are  described.  The  model  velocity  ranges  between 
4,  000  and  17,  000  ft/sec  with  range  pressures  varying  from  15  to  380  mmHg 
in  air.  The  experimental  data  is  discussed  in  detail,  and  the  analysis  is  used 
to  compute  and  correlate  appropriate  transition  parameters. 

Based  on  the  correlation,  predictions  of  the  location  of  transition  in  the  wake 
of  a  12-degree  half-angle  cone  at  22,  000  ft/ sec  are  made  for  the  altitude  range 
from  200,  000  to  100,  000  feet.  For  example,  transition  is  placed  at  approximate 
ly  250,  100,  and  20  feet  from  the  apex  of  a  5 -foot  cone  at  altitudes  of  200,  000, 
150,  000  and  100,  000  feet,  respectively,  assuming  the  boundary  layer  on  the 
cone  surface  remains  laminar. 
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NOMENCLATURE 


Symbols 

Fpl 

h 

H 

j 

I 

M 

N 

P 


Pr 


r.  y 


R 

Re 

X 


Rb-  d 


s 

T 

u 

V 

V 

X 


p 


Integrals  defined  by  equations  (6)  and  (7) 
static  enthalpy 
stagnation  enthalpy 
an  index,  either  0  or  1 

PP/PI  Ml 
Mach  number 

number  of  strips 

pressure 

Prandtl  number 

normal  coordinate,  measured  from  axis  or  plane  of  symmetry 
gas  constant 

Reynolds  number,  based  on  local  inviscid  properties  and 
X  distance. 

base  radius  of  body  or  diameter,  respectively 

Sutherland  constant 

static  temperature 

x-component  of  velocity 

r-component  of  velocity 

velocity  vector 

streamwise  coordinate,  measured  from  forward  stagnation 
point 

mass  density 
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NOMENCLATURE  (Concl'd) 


Symbols 

( 

V 

s 

e 

X 

y 

r 

o 

Subscripta 

1 

c 

k 

*.  t.  »? 


viacoaity 

tranaformed  atreamwiae  coordinate 

tranaformed  normal  coordinate 

wake  thickneaa  (from  axia  or  plane  of  aymmetry) 

momentum  thickness,  or  cone  angle 

transformed  thickness  parameter 

ratio  of  specific  heats 

transformed  shear  function 

transformed  heat  diffusion  function 

evaluated  at  r  =  S  ox  r)  -  1 
evaluated  at  r  =  o  or  77  =  0 
evaluated  at  edge  of  kth  strip 

partial  derivative  with  respect  to  variable  indicated 
reservoir  conditions 
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I.  INTRODUCTION 


The  flow  in  the  wake  of  hypersonic  vehicles  has  been  the  subject  of  consider¬ 
able  research  in  recent  years.  The  interest  in  this  subject  stems  mainly  from 
the  realization  that  the  trails  of  hypersonic  vehicles  are  a  potentially  signifi¬ 
cant  source  of  observables,  i.  e. ,  electrons  and  radiating  species.  The  bulk 
of  the  lieterature  deals  mainly  with  the  trails  of  blunt  bodies  2,  3  where  the 
viscous  core  comprises  a  significant  portion  of  the  wake  flowfield  (from  the 
point  of  view  of  observables)  only  if  turbulent  mixing  is  present.  However, 
in  the  case  of  slender  bodies,  the  observables  in  the  inviscid  flowfield  are 
minimized  and  the  viscous  region  (laminar  or  turbulent)  becomes  the  sole 
source  of  electrons  and  radiating  species. 

An  excellent  account  of  the  state  of  the  art  of  the  turbulent  far  wake  has  been 
presented  by  Lees^  and  Vaglio-Laurin,  Bloom,  and  Byrne.  ^  These  authors 
extended  their  blunt-body  analyses  to  the  slender-body  case.  Most  analyses 
of  the  laminar  far  wake  (shown  schematically  in  figure  1)  have  been  restricted 
to  similar  solutions  involving  linearized  forms  or  to  solutions  obtained  by  the 
Karman-Pohlhausen  method.  The  notable  exception  is  the  classical  Goldstein 
solution  for  the  incompressible  wake  of  a  flat  plate,  ^ 

Recent  experimental  data  acquired  at  Avco  RAD  and  the  Naval  Ordinance 
Laboratory  (NOL)  ballistic  ranges  show  appreciable  lengths  of  laminar  wake 
at  Reynolds  numbers  higher  than  those  previously  reported.^  Therefore,  for 
slender  bodies  a  critical  analysis  of  the  laminar  far  wake  has  a  twofold  value: 
(1)  to  provide  detailed  flowfields  for  the  prediction  of  observables,  and  (2)  to 
compute  accurately  the  appropriate  transition  parameters. 

This  paper  presents  a  theoretical  analysis  of  the  laminar  far  wake  and  also 
contains  some  of  the  experimental  results  obtained  in  the  Avco  RAD  ballistic 
range.  The  method  of  theoretical  analysis  is  an  extension  of  the  multistrip 
solution  of  the  laminar  boundary-layer  equations.^  No  assumptions  about  the 
form  of  the  solution  are  required;  only  the  initial  and  boundary  conditions  must 
be  specified.  Since  no  solution  of  the  near  wake  is  presently  available,  the 
initial  profiles  are  somewhat  arbitrary.  However,  initial  velocity  and  enthal¬ 
py  profiles  can  be  prescribed  with  some  degree  of  confidence  from  existing 
knowledge  of  gross  near -wake  properties  and  by  utilizing  overall  momentum 
and  energy  conservation.  Only  pure  air  wakes  in  chemical  and  thermodynamic 
equilibrium,  or  completely  frozen,  are  considered  in  the  present  analysis. 
Some  results  for  nonequilibrium  wakes  have  been  obtained  by  Bloom,  et  al. , 
by  the  Karman-Pohlhausen  technique.  ^  The  nonsimilar  analysis  with 
chemical  nonequilibrium  has  been  formulated  and  some  preliminary  resvilts 
obtained  but  they  are  not  included  here. 
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II.  NONSIMILAR  ANALYSIS 


Steady  flow  downstream  of  the  neck  of  a  wake  (i.  e.  ,  the  stagnation  point  in  the 
wake  beyond  which  no  recirculation  takes  place)  can  be  described  in  the  context 
of  boundary-layer  theory  by: 


(pufi)^  +  (pvr))j  =  0 

1  i 

puUj  +  pvUj  =  -  Pjj  +  —  (pr'Uj)j 


where 


(1) 

(2a) 

(2b) 

(3) 


0  for  planar  flow 


1  for  axisyffltnetric  flow. 


The  conventional  boundary  conditions  apply  at  the  outer  edge  of  the  wake, 


8u 

u  «  u  1  and  — —  =  0 
^  dr 


dH 

H  -  H,  and  -  =  0 

‘  dr 


I  at  r  = 


at  r  =  5  , 


while  at  the  plane  or  axis  of  symmetry  the  boundary  conditions  are 


du 

u  »  u-  and  —  =  0 

dr 


dH 

H  -  H-  and  — —  -  0 

‘  dr 


at  r  9:  0 


Parenthetically!  it  may  be  noted  at  thie  point  that  if  P,  >  0  and  Pr  ■  1,  the 
energy  equation  it  catiefied  by  the  Crocco  integral; 

H  -  Au  B 


where  a  and  b  are  conetanta: 


is  a  particular  solution  of  the  energy  equation.  This  solution  is,  of  course, 
only  valid  if  the  initial  conditions  also  satisfy  the  Crocco  integral.  The  present 
analysis  is -not  restricted  to  this  case,  since  the  Crocco  integral  is  not  em¬ 
ployed. 

In  the  present  analysis,  the  wahe  is  divided  into  an  arbitrary  number  of  strips, 
N,  in  the  streamwise  direction,  and  the  equations  then  are  integrated  with 
respect  to  r  across  the  strips,  i.  e. ,  between  the  limits  of  t  =  0  and  r  =  r^, 
to  give  two  ordinary  differential  equations:* 


*Th«  itrip-iattsnl  method  ie  diecuesed  in  gteeter  detail  in  reference  7. 
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The  transformed  independent  variables  rj,  are  now  introduced 


where  L  is  a  characteristic  length  of  the  problem  and 


1 


is  the  thickness  of  the  velocity  and  thermal  profiles  in  the  transformed  plane. 
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The  equations  now  become 
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I 


PM 
Pi  Ml 


The  transformed  boundary  conditions  are 
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Moreover!  the  transformed  momentum  and  energy  equations  evaluated  at 
>1*0  become'" 


The  velocity  and  stagnation  enthalpy  profiles  can  be  expressed  as  polynominals 
in  7.  However!  additional  numerical  accuracy  can  be  obtained  by  choosing  the 
following: 


N 


n  •  0 


N+  I 


which  identically  satisfy  the  imposed  boundary  conditions  at  »?  =  1  (and  moreover 
satisfy  the  requirement  that  all  higher-order  derivatives  also  vanish  at  that 


*Noct  that 


t-*0  ,  ^"t»0 


8- 


point).  The  polynominal  coefficients  are  made  to  satisfy  the  boundary  condi¬ 
tions  at  i;  B  0  and  the  values  of  u](  and  H](  at  the  strip  boundaries  and  at  the 
wake  axis,  or  plane  of  symmetry,  and  also  Hrjric  ■ 

The  ordinary  equation  of  state.  P  ■  pRT.  provides  the  necessary  thermodynamic 
relations  for  a  perfect  gas.  while  curvefits  of  the  Mollier  chart^^  can  be  em¬ 
ployed  for  equilibrium  air. 

Frozen  flow  can  be  approximated  by  the  use  of  the  perfect  gas  law  with  appro¬ 
priate  gas  constant.  R.  and  ratio  of  specific  heats,  y.  At  low  temperatures, 
the  Sutherland  viscosity  law  is  employed  to  give 


where 

S  =  120°K 


while  at  elevated  temperatures  the  pp  ratio  is  approximated  by 


where 


A  =  3.  03 
B  =2.03 

hj  =  4.  5  X  10^  ft^/  sec2. 


Therefore.  2N  +  2  first-order  ordinary  differential  equations  are  obtained 
(N  momentum  integral.  N  energy  integral,  and  momentum  and  energy  at  t  =  0) 
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which  can  be  eolved  numerically  for  the  2N  +  2  unknowns  (A,  for  k  «  2, 

3» . . . ,  N,  and  I  and  ),  A  predictor-corrector  marching  scheme  is 
employed  for  solution  of  the  equations  on  an  IBM  7094.  * 


*It  is  impoccant  to  note  that  in  a  numericsi  solution  equations  (8)  and  (9)  ate  singular  for  =  0,  although  as  pointed  out 

by  T,  Kubota  (in  pciTatc  communication),  this  difficulty  can  be  avoided  in  an  analytical  solution  by  use  of  an  appropti- 
ate,  nonseto  pressure  gradient.  The  behavior  of  the  equations  nest  this  tear  stagnation  point  has  been  examined  oy 

S.l.  Cheng  (Avco  RAD-TM^3<23)  and  by  Vaglio-Lautin,  Bloom  and  Byrne^.  In  the  present  analysis,  solution  of  the  equa¬ 
tions  is  always  begun  slightly  away  from  the  stagnation  point,  with  a  small  but  finite  value  of  u^. 
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m.  COMPARISON  OF  METHODS  OF  APPROACH  TO  THE  FLUID 
DYNAMICS  OF  THE  WAKE 


The  laminar  wake  of  a  slender  body  only  recently  has  gained  practical 
importance  as  a  fluid-dynamic  problem  due  to  its  apparent  high  stability  at 
hypersonic  speeds’**  which  has  not  been  previously  observed,  Howeveri  compu¬ 
tations  of  the  laminar  wake  of  a  thin  flat  plate  at  low  speeds  were  reported  in 
the  1930' s  (e.  g.  •  reference  6).  and  various  approximate  methods  have  been 
developed.  The  only  "exact"  solution,  however,  is  that  of  Goldstein,  ^  which  is 
essentially  a  finite  difference  solution  for  the  wake  in  the  immediate  vicinity 
of  the  trailing  edge.  The  approximate  methods  which  have  been  developed  fall 
into  essentially  two  categories-- (U  the  "similarity"  solutions,  and  (2)  the 
Karman-Pohlhausen-type  solutions. 

The  similarity  analyses  require  linearization  of  the  velocity  defect  (uj  .  u) 
and  usually  neglect  the  streamwise  pressure  gradients,  which  results  in  a 
closed-form  solution  of  the  partial  differential  equations  of  momentum  and 
energy  satisfying  the  wake -boundary  conditions.  The  effect  of  initial  conditions 
is  felt  only  through  two  constants  in  the  solution,  which  are  usually  defined 
by  overall  momentum  and  energy  conservation  from  the  body  boundary  layer. 

The  Karman-Pohlhausen-type  analysis,  as  the  name  implies,  is  an  extension 
of  the  integral  boundary -layer  approach  to  the  wake  problem.  Clearly,  the 
accuracy  of  an  integral-type  solution  is  dependent  on  its  ability  to  represent 
the  detailed  behavior  of  properties  across  the  layer,  which  is  a  function  of 
the  number  of  independent,  free  parameters  retained  in  the  solution.  **  In 
the  Karman-Pohlhausen  method,  the  only  free  (undetermined)  parameters  are 
usually  the  thickness  S,  the  velocity,  and  enthalpy  on  the  axis,  or  plane,  of 
symmetry,  uc  and  ,  and  the  curvature  of  the  enthalpy  profile  at  the  axis 

.  Hence,  the  only  nonsimilar  influence  on  the  normalized  velocity  profile 

9  the  thickness  3,  and  the  only  influence  on  the  normalized  enthalpy 

- )  is  due  to  5  and  the  curvature  of  the  profile  at  the  axis.  The 

Hj-Hcy 

initial  conditions  are  specified  through  the  initial  values  of  the  free  parameters. 
The  momentum -integral  equation  itself  is  no  more  than  an  expression  of  the 
fact  that  the  drag  contained  in  the  wake  is  constant,  implying  that  in  a  numerical 
integration  of  the  equations  the  initial  conditions  must  reflect  the  correct  drag 
(i.  e.,  initial  momentum  thickness)  to  guarantee  the  proper  integration  constant. 
In  the  Karman-Pohlhausen  technique,  this  requirement  can  cuase  a  mismatch 
of  the  overall  thickness,  depending  on  the  shape  of  the  initial  profile  being 
approximated. 


"rV 

profile  I 


'Ditciwtid  in  futthet  detail  later. 

* ‘Reference  8  is  an  excellent  example  of  the  use  of  the  Karman-Pohlhausen  technique  in  studying  nonequiUbrium  viscous 
flow  problems. 
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The  present  method  aleo  ie  categorized  as  an  integral  method  (being,  in  fact, 
an  application  of  Oorodnitsyn'e  "method  of  integral  relatione")  and.  ae  such, 
is  an  approximate  method  also.  However,  it  avoids  the  restrictions  associated 
with  the  Karman-Pohlhausen  method  since  the  velocity  and  enthalpy  on  the 
strip  boundaries  are  additional  free  parameters  which  must  be  determined 
from  the  differential  equations. 


illustrated  in  figures  2  and  3.  The  similarity  solutions  give 


The  problem  of  the  laminar  wake  of  a  thin  flat  plate  was  examined  in  detail 
to  compare  the  approximate  methods  mentioned  above  and  the  present  method  with 
the  classical  Goldstein  solution.  This  problem  also  has  the  feature  that  there 
is  no  recirculation  and  the  initial  (Blasius)  profile  is  known  exactly.  The  decay 
of  the  velocity  along  the  plane  of  symmetry  and  typical  velocity  profiles  are 

*  ..  ..  .  ..  .  .^l/jl/2 

which  is  particularly  poor  near  the  trailing  edge,  although  the  rate  is  correct 
asymptotically  at  large  distances.  The  Karman-Pohlhausen  method  shows 
remarkably  good  agreement  with  the  exact  results  near  the  center  of  the  wake 
in  this  Cwse.  Some  differences  in  thickness  S  are  illustrated  in  figure  4. 

The  present  method,  however,  is  in  essentially  perfect  agreement  with 
Goldstein's  solution.  It  should  be  noted  that  the  results  from  the  present 
method  shown  were  obtained  with  six  strips,  but  calculations  also  were  performed 
with  two  through  five  strips.  Rapid  convergence  of  the  solution  was  obtained 
with  greater  than  four  strips,  although  two  strips  offered  little  improvement 
over  the  Karman-Pohlhausen  (one  strip)  solution. 


'The  tefiieoee  length  L  ie  taken  ne  the  plate  length  in  thie  case. 
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Figure  3  VELOCITY  PROFILES  IN  THE  WAKE  OF  A  FLAT  PLATE 

63-3488 


L 


Figur*  4  GROWTH  OF  THE  WAKE  OF  A  FLAT  PLATE 
63-4980 


IV.  HYPERSONIC  LAMINAR  FAR  WAKE  FOR  A  SLENDER  BODY 


A.  INITIAL  (NECK)  CONDITIONS 

The  conditions  which  occur  at  the  neck  of  the  hypersonic  wake  are  the  result 
of  simultaneous  mixing  and  recirculation  of  the  boundary  layer  shed  from  the 
body,  as  well  as  a  strong  adverse  pressure  gradient  (which  "drives"  the  re¬ 
circulating  flow).  Also,  the  base  of  the  body  is  a  stagnation  point  which  has 
been  experimentally  observed  to  have  many  of  the  features  (e.  g.  ,  high  heat 
flux)  of  a  blunt-body  stagnation  point,  although  less  severe.  Unfortunately, 
the  near  wake  has  thus  far  eluded  successful  theoretical  treatment,  but  it  is 
currently  receiving  considerable  attention.  Various  approximations  have  been 
attempted,  but  none  offer  any  really  satisfactory  insight  into  the  mechanism 
or  magnitude  of  the  recirculation  and  base  heating. 

The  mixing  phenomenon  per  se  has  been  treated  by  Denison  and  Baum,  ^^Vaglio- 
Laurin,  Bloom  and  Byrne  ^  as  an  extension  of  the  earlier  ideas  of  Chapman, 
who  studied  the  mixing  of  a  stream  with  quiescent  air.  Chapman's  analysis 
assumed  a  zero-thickness  shear  layer  at  the  beginning  of  the  mixing  region, 

u* 

for  which  he  obtained  a  value  of  —  =  0.  587  for  the  velocity  ratio  of  the 

“1 

dividing  streamline.  Denison  and  Baum  considered  the  supersonic  mixing 
process  starting  with  a  Blasius  profile  and  both  cold  quiescent  air  (free -stream 
static  temperature)  and  hot  quiescent  air  (free -stream  stagnation  temperature) 
for  a  variety  of  cone  and  wedge  angles.  Their  results  indicate  that  the  dividing 
streamline  velocity  ratio  is  on  the  order  of  0.  25  to  0.  20  at  hypersonic  speeds 
for  10-to  20-degree  cones  for  the  cold  case,  and  0.  35  to  0.  30  for  the  hot  case. 
With  constant-temperature  quiescent  air,*’"’‘con8tant  pressure  mixing,  and 
Prandtl  number  of  unity,  the  Crocco  integral  is  a  solution  of  the  energy  equa¬ 
tion,  giving: 


where 

H*  =  dividing  streamline  total  enthalpy 
H(j  =  total  enthalpy  of  the  quiescent  air 

=  total  enthalpy  of  the  inviscid  stream. 

*The  itreamlinc  which  tepkratei  the  flow  which  peases  the  neck  from  the  recirculating  flow. 
••At  apptoximacely  the  body  wall  temperature. 
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Hv  H*  •  *** 

For  th«  cold  cas«i  I.  «  li  giving  —  «  —  •  Clearly,  in  the  hot  case  |r  ■  i* 

Hi  Hi  “1  ^ 

Therefore,  coneideration  of  the  mixing  process  alone  leads  to  a  realistic  lower 
limit  on  tlie  dividing  streamline  enthalpy.  Assuming  that  the  compression  pro¬ 
cess  occurs  through  a  shocklike  discontinuity  (allowing  the  neglect  of  viscous 
dissipation),  then  the  dividing  streamline  enthalpy  can  be  taken  as  the  value 
occurring  on  the  axis  at  the  neck  (i.  e.  ,  where  the  dividing  streamline  is  stag¬ 
nated). 

The  nature  of  the  compression  process  in  the  neck  region  suggests  that  the  neck 
bears  certain  similarities  to  a  boundary -layer  separation  point.  (Indeed,  it 
has  been  recognized  by  various  authors  that  gross  properties  such  as  base 
pressure  and  wake  angle  are  predictable  in  terms  of  separation-point  pressure 
correlations.  )  Therefore,  the  use  of  a  separation -point  velocity  profile  from 
similarity  solutions  such  as  Cohen  and  Reshotko^^  appears  to  be  a  reasonable 
approximation.  However,  the  enthalpy  profile  which  would  result  after  such  a 
compression  is  not  known  beyond  the  axis  value  indicated  above.  Clearly,  the 
use  of  the  Crocco  integral  is  not  correct  for  the  entire  near  wake,  a  region 
which  Includes  strong  pressure  gradients;  but,  it  does  have  the  justification  of 
being  at  least  a  consistent  means  of  approximating  the  neck  enthalpy  profile  in 
terms  of  the  assumed  velocity  profile  and  the  (approximately)  known  enthalpy 
on  the  axis.  However,  it  should  be  noted  that  the  use  of  the  Crocco  integral 
leads  to  an  "annulus"  of  high  temperature  gas  rather  than  a  "cylinder"  at  the 
neck  for  slender  bodies.  It  is  not  obvious  whether  this  is  physically  reason¬ 
able;  intuitively,  it  would  be  expected  that  the  peak  temperature  would  occur 
at  the  axis.  Regardless,  for  the  lack  of  any  better  knowledge  of  the  near  wake, 
the  Crocco  integral  will  be  used  in  the  present  analysis  to  specify  the  initial 
total  enthalpy. 

The  necessity  of  matching  the  initial  drag  (i.  e.  ,  momentum  thickness)  at  the 
neck  has  been  previously  mentioned.  Assuming,  a  priori,  that  the  drag  added 
in  the  near  wake  is  negligible  in  comparison  with  that  shed  from  the  body  in  the 
boundary  layer,  the  initial  momentum  thickness  is  given  by: 


[pi  ui^  "  tPl“l^0+l)R^l9)body 


where 


^neclt  •  <(> 


1/i  + 1 
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In  the  numerical  examples  the  initial  value  of  the  thickness  parameter,  X 
computed  from  the  neck  momentum  thickness  in  the  following  manner: 


A 


I  n  . F,_  ^ 


2/j  +1 


is 


The  above  assumptions  are  verified  a  posteriori  only  in  that  they  lead  to  results 
which  are  in  good  agreement  with  the  experimentally  observed  facts  (i.  e.  ,  the 
wake  thickness  and  its  Reynolds -number  dependence),  which  is  discussed  in 
greater  detail  later  in  this  paper.  However,  when  they  are  applied  at  very  low 
pressures  (high  altitude)  certain  anomalous  results  are  obtained,  namely  initial 
wake  thicknesses  which  exceed  the  body  diameter  and  even  the  "displaced" 
body  diameter  (i.  e,  ,  body  plus  boundary  layer).  This  is  undoubtedly  due  to 
neglect  of  interaction  with  the  individual  flow.  The  viscous -inviscid  interaction 
is  a  separate  problem  which  is  not  considered  in  this  analysis,  nor  in  any  of  the 
examples. 

B.  NUMERICAL  EXAMPLES 

The  present  analysis  was  applied  to  study  the  development  of  the  laminar  far- 
wake  for  some  typical  cases  of  interest.  The  medium  is  assumed  to  be  air  in 
equilibrium,  unless  otherwise  noted,  with  a  constant  Prandtl  number  of  0,72, 

The  following  case  is  an  example  which  is  typical  of  experiments  performed  in 
the  Avco  RAD  ballistics  range: 

10-degree  (semivertex  angle)  cone 
base  diameter:  0.  300  inch 
velocity:  15,  688  ft/sec 
pressure:  150  mmHg 

He 

Velocity  and  temperature  profiles  (with  —  -  0.  3  initially)  for  the  wake  at 

these  conditions  are  shown  in  figures  5  and  6,  The  features  of  the  flowfield 
which  can  be  observed  from  the  Schlieren  photos  (e,  g,  ,  figure  14,  the  thickness 
and  growth  rate)  are  essentially  in  agreement  with  these  results. 

This  slender  axisymmetric  body  example  may  be  contrasted  to  the  results  for 
a  blunt  semi-infinite  body;  namely,  a  cylinder  placed  normal  to  the  flow.  The 
conditions  selected  correspond  to  one  of  the  experiments  performed  by  McCarthy 
in  the  GALCIT  with  tunnel: 

0.  300 -inch-diameter  cylinder 
free-stream  Mach  number:  5.  7 
reservoir  pressure:  35  psig 

reservoir  temperature:  262®F 
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Figure  5  VELOCITY  PROFILES  IN  THE  EQUILIBRIUM  WAKE  OF  A  10-DEGREE  CONE 

63-3501 


In  this  cassi  the  viscous  core  may  entrain  mass  flow  from  the  entropy  layer  as 
it  grows.  This  "swallowing"  effect  was  estimated  using  the  method  of  character¬ 
istics.  It  was  found,  however,  that  the  laminar  core  (in  the  present  example) 
does  not  entrain  any  appreciable  mass  from  the  entropy  layer  for  at  least  25 
diameters  from  the  neck.  An  appreciable  velocity  gradient  is  encountered, 
nonetheless,  due  to  the  presence  of  expansion  waves  reflected  from  the  bow 
shock. 

A  constant  stagnation  enthalpy  (equal  to  the  free -stream  value)  was  assumed, 
since  the  tunnel  running  time  was  long  allowing  the  model  to  approach  an  equili¬ 
brium  temperature.  Typical  resulting  velocity  and  temperature  profiles  are 
shown  in  figures  7  and  8.  The  very  fast  velocity  decay  observed  here  is  similar 
to  that  found  in  the  classical  flat-plate  solution,  and  appears  to  be  peculiar  to 
semiinfinite  bodies.  The  relative  importance  of  the  nonsimilar  behavior  of 
the  solution  and  also  of  the  attendant  favorable  pressure  gradient  in  this  ex¬ 
ample  are  illustrated  in  figure  9.  As  much  as  a  70 -percent  error  can  be  attribut¬ 
ed  to  neglect  of  the  pressure  gradient.  *  It  should  also  be  noted  that  no  effect 
of  Reynolds  number  on  the  decay  rates  can  be  expected  (also,  as  in  the  flat- 
plate  case)  if  the  initial  value  of  the  thickness  parameter  \  is  constant,  as  indi¬ 
cated  by  McCarthy's  data. 

As  an  example  for  the  laminar  far  wake  of  a  slender  body  at  flight  conditions, 
a  constant  velocity  (22,  000  ft/sec)  trajectory  from  200,  000  to  100,  000  feet  was 
assumed  for  a  10 -foot,  12-degree  (semivertex)  angle  cone.  An  initial  enthalpy 
H 

ratio  of  —  =  0.  3  was  assumed  and  the  calculations  were  performed  for  both 

«1 

frozen  (  y  =1.4  and  R  =  1,  724  ft^/sec^  /  *  r  )  and  equilibrium  air.  Typical 
profiles  and  the  velocity  and  temperature  along  the  wake  axis  are  given  in 
figures  10,  11,  12,  and  13.  The  results  are  qualitatively  similar  to  those  for 
the  previous  slender -body  example.  It  is  interesting  to  note  the  growth  in 
length  of  the  trail  with  decreasing  altitude.  However,  at  200,  000  feet  the  ini¬ 
tial  thickness  exceeds  the  body  diameter  by  about  30  percent.  Consideration 
of  the  interaction  phenomena  should  lead  to  a  shorter  trail  at  this  altitude. 

(The  equilibrium  results  at  200,  000  feet  and  the  frozen  results  at  100,  000  feet 
are  included  only  for  comparison's  sake;  clearly  neither  are  physically  realistic.  ) 


'Attention  muit  be  called  to  the  feet  that  much  stnaller  preaauie  gradients  are  encountered  with  slender,  sharp  cones. 
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Figure  9  VELOCITY  DECAY  IN  THE  WAKE  OF  A  CYLINDER 
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Figure  10  VELOCITY  DECAY  ALONG  THE  WAKE  AXIS  OF  A  lO-FOOT,  12-DEGREE  CONE 


-28 


Figure  12  VELOCITY  PROFILES  IN  THE  WAKE  OF  A  lO-FOOT,  12-DEGREE  CONE 


Figure  13  TEMPERATURE  PROFILES  IN  THE  WAKE  OF  A  lO-FOOT,  12-DEGREE  CONE 


V.  EXPERIMENTAL  STUDIES 


A  series  of  ballistic -range  experiments  with  conical  models  of  10,  15,  and 
27-1/2-degree  half  angles  were  performed.  The  model  velocities  range  be¬ 
tween  4,000  and  17,000  ft/sec  with  range  pressures  varying  from  15  to 
380  mmHg  in  air. 

The  techniques  for  the  use  of  a  light-gas  gun  for  launching  projectiles  to 
reentry  velocities  are  well  established.  The  projectile  must  have  aerodynamic 
stability  and  must  be  capable  of  withstanding  the  launch  loading  of  the  gun.  A 
conical  model  with  a  semivertex  angle  between  10  and  20  degrees  has  been  de¬ 
veloped  which  can  withstand  the  launching  environment  of  a  light -gas  gun  to 
achieve  velocities  up  to  17,000  ft/sec.  This  model,  shown  in  figure  14,  consists 
of  a  steel  tip  and  titanium  afterbody.  It  has  a  hollow  base  to  provide  a  static 
margin  of  15  percent.  The  sabot  consists  of  an  aluminum -driven  disc  and 
four  Zel  IX  petal  components  that  interlock  with  the  model  to  provide  the  neces¬ 
sary  inbore  support  of  the  model  during  launch.  Many  sabot-model  combinations 
were  tried  before  this  successful  package  was  achieved.  Failure  of  the  model 
was  indicated  by  fracture  and  separation  of  the  tip  or  by  bulging  or  separation 
at  the  bimetal  interface.  These  models  have  a  dispersion  at  the  end  of  the 
range  (approximately  40  feet  from  the  muzzle)  of  about  0.  75  inch.  Approxi¬ 
mately  90  percent  of  these  models  result  in  successful  launchings.  The  pitch 
and  yaw  of  the  models,  as  measured  from  the  photographic  stations  of  the 
range,  are  generally  limited  to  a  few  degrees  as  measured  from  the  velocity 
vector. 

A  one-piece,  steel-cone  model  (having  a  hollow  base)  has  been  used  for  half 
angles  of  27.  5  degrees.  These  models  exhibit  a  larger  dispersion  than  the 
above  two-piece  design. 

Two  0.  600 -caliber  light-gas -gun  ranges  have  been  used  in  this  program. 

The  first  range  is  equipped  with  a  12-inch  diameter,  double -traverse  Schlieren 
system.  The  12-inch  diameter  parabolic  mirror  is  located  inside  a  10-foot- 
diameter  by  14-foot-long  tank.  This  range  has  been  used  to  study  the  low- 
velocity  regime  (under  11,  000  ft/sec)  where  the  sensitivity  of  the  double¬ 
traverse  system  is  utilized  for  flow  visualization  at  low  density  (e.  g.  ,  pres¬ 
sures  in  the  range  15<  P<100  mmHg),  A  0.0005-inch  Mylar  foil  screen  has 
been  used  to  provide  the  trigger  pulse  for  the  spark  light  source  in  the  Schlieren 
station  in  some  of  the  photographs.  This  screen  was  installed  to  obtain  greater 
reliability  over  the  light  screen  originally  employed  within  the  vacuum  tank  and 
to  permit  easy  adjustment  of  the  position  of  the  model  along  its  trajectory  when 
the  photograph  is  made.  There  is  no  effect  from  the  foil  on  the  flowfield  since 
there  is  negligible  time  for  any  disturbance  to  propagate  before  the  photograph 
is  made. 

•A  mote  detailed  description  of  the  Avco  RAD  ballistics-range  facility  and  these  tests  is  contained  in  Avco  RAD-TM-63-20. 
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A  second  caliber -0. 600  light-gas-gun  range  has  been  used  to  obtain  data  from 
11,  000  to  17,  000  ft/sec  at  pressures  greater  than  30  mmHg  in  air.  Primary 
instrumentation  on  this  range  is  a  12-inch  diameter,  single-pass  Schlieren 
apparatus. 

Some  results  covering  the  range  of  experimental  conditions  for  which  transi¬ 
tion  was  clearly  observable  from  the  Schlieren  photographs  are  compiled  in 
table  1.  Five  of  the  Schlieren  photographs  which  exemplify  the  typical  features 
of  the  flowfield  of  a  slender-body  wake  are  presented  here  as  figures  15,  16, 

17,  18,  and  19.  Similar  features  can  be  seen  in  shadowgraph  photos  (e.  g.  , 
figure  1,  reference  5).  It  should  be  noted  that  the  near -wake  details  cannot  be 
seen  in  such  photographs.  Rather,  recourse  must  be  made  to  observables  such 
as  the  apparent  origin  of  the  trailing  shock  and/or  a  minimum -thickness  section 
to  define  the  neck  location. 

Model  pitching  in  the  plane  of  the  photograph  is  clearly  identifiable  by  the  asym¬ 
metric  location  of  the  wake  within  the  shock  layer.  Small  yaw  angles  normal  to 
the  plane  of  the  photograph  are  indicated  by  an  abnormally  thick  wake  as  com¬ 
pared  to  zero  yaw  at  the  same  conditions,  and  by  appearance  of  a  streak  down 
the  center  of  the  wake. 

The  occurrence  of  turbulence  is  always  clearly  identifiable  for  sharp  bodies* 
and  is  usually  preceded  by  an  instability  which  manifests  itself  as  sinusoidal 
waviness  of  the  wake.  Generally,  the  waviness  is  followed  almost  immediately 
followed  by  turbulence,  but  occasionally  it  persists  for  an  appreciable  length 
before  breakdown  into  turbulence.  In  none  of  the  photos  are  the  waves  observed 
to  be  opposing  (i.  e. ,  peak  opposite  peak,  valley  opposite  valley).  Some  evidence 
of  helical  motion  can  be  seen  in  the  cases  with  yaw  (e.  g.  ,  figure  15).  By  con¬ 
trast,  a  model  without  yaw  at  approximately  the  same  conditions  is  shown  in 
figure  16.  The  effect  of  pressure  on  transition  is  demonstrated  in  figures  17, 

18,  and  19. 


•This  is  not  necessarily  the  esse  for  blunt  bodies,  where  the  viscous  core  may  be  obscured  by  the  inviscid  entropy  layer 
which  enshrouds  it. 
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Figure  15  TYPICAL  FEATURES  OF  A  SLENDER  BODY  FLOW  FIELD.  EXAMPLE  1 
10*D«gr««  Cona  V  =  14.800  ft>sec 

0.300-Inch  Bosa  Diamatar  P  =  150  mm  Hg 
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Figure  18  TYPICAL  FEATURES  OF  A  SLENDER  BODY  FLOW  FIELD,  EXAMPLE  4 
15'Degree  Cone  V  =  10,000  ft-sec 

0.400-lnch  Bose  Diameter  P  =  250  mm  Hg 
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AL  FEATURES  OF  A  SLENDER  BODY  FLOW  FIELD,  EXAMPLE  5 
)-Degr*a  Con*  V  =  9,000  ft-sec 

,300>1nch  Bom  Diomater  P  =  350  mm  Hg 


VI.  TRANSITION  CORRELATION 


It  is  generally  agreed  that  wake  transition  is  dependent  on  a  Reynolds  number 
and  a  Mach  number  which  are  based  on  the  relative  velocity  of  the  wake  (uj  - 
Uc  ),  (references  4,  19.  20,  and  21).  Stability  theory  indicates  that  the  Reynolds 
number  and  Mach  number  should  be  evaluated  using  physical  properties  at  the 
"critical  point"  in  the  layer.  *  However,  transition  correlations  commonly  use 
properties  at  the  outer  edge  of  the  layer  and  introduce  a  temperature  ratio  as  a 
third  parameter.  Several  length  scales,  namely  surface  distance  to  transition 
point,  boimdary-layer  thickness,  momentum  thickness,  etc.  ,  have  been  employed 
in  defining  a  transition  Reynolds  number.  None  of  these  parameters  are  com¬ 
pletely  satisfactory  in  situations  where  the  viscous -layer  history  is  significant 
(e.  g.  ,  strongly  varying  edge  properties,  nonsimilar  profile  shape,  etc.). 


Clearly,  the  latter  remarks  are  appropriate  to  wake  transition.  A  heuristic 
approach  to  this  problem  (short  of  a  complete  stability  analysis)  is  to  attempt 
correlation  using  physically  reasonable  length  from  the  origin  of  the  viscous 
layer  to  the  transition  point,  ,  and  the  layer  thickness  at  transition, 
both  of  which  are  readily  measurable.  A  third  possibility,  however,  is  the 
relative  distance*  *=I*a  particle  travels  along  the  axis  to  the  transition  point,  . 
This  is  a  measure  of  the  time  which  a  disturbance  which  is  convected  with  the 
particle  has  to  amplify.  The  significance  of  this  time  is  suggested  by  the 
experiments  of  Sato  and  Kuriki,^^  in  which  the  disturbances  were  observed  to 
originate  at  the  center  of  the  wake  and  grow  as  they  convected  downstream. 
The  distance  is  given  by: 


X 


(uj  -  u^)  dt 


(reference  system  fixed  with  respect  to  the 
vehicle). 


A  photograph  of  the  flowfield,  however,  only  reveals  the  relative  distance  travel¬ 
ed  by  a  particle  along  the  axis  as  compared  with  the  distance  traveled  by  the 
vehicle: 


^^tr^photo  “  ^ 


Up  dt 


*  Gold,  H.,  Avco  RAD,  Personal  Commuaication 
is  defined  as  the  woke  diameter  at  transition. 

•••i.e.,  in  a  coordinate  system  fixed  with  respect  to  the  medium  in  which  the  body  is  immersed. 
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Therefore, 


dx 


where 


dt 


and 


^^tr^photo 


L 


The  distance  to  transition  along  a  wall  or  in  a  shear  layer  is  consistent  with 
this  formulation,  (i.  e.  ,  with  “c  =  0)* 

The  preceeding  analysis  was  used  to  correlate  the  wake -transition  data  from 
the  experiments  in  the  Avco  RAD  ballistics  range  and  some  from  the  NOL 
ballistics  range.  The  Reynolds  number  and  Mach  number  are  defined  as: 


The  values  of  —  and  Xj,  were  computed  for  the  experimental  conditions 
“l 

utilizing  the  values  of  measured  from  Schlieren  and  shadowgraph  photos. 

The  value  of  5(,  was  also  measured.  (Transition  was  defined  by  the  occurrence 
of  waviness  at  the  outer  edge  of  the  wake.)  A  correlation  of  these  experimental 
data  is  offered  in  figures  20,  21,  and  22,  based  on  X^  ,  8^^  ,  and  X^^  ,  re¬ 

spectively.  The  scatter  is  at  least  partially  due  to  the  fact  that  no  attempt  to 
accovmt  for  temperature  effects  has  been  included,  and  the  ratio  T^/Tj  at  tran¬ 
sition  covers  a  range  from  about  4  to  10  for  the  slender,  high-velocity  cones  based 
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He 

on  an  assumed  value  of  ■“  =  0.  3  at  the  neck. 

Hi 

15  1 

Some  transition  data  of  Sato  and  Kuriki  and  McCarthy^®  also  are  included 
for  the  sake  of  comparison.  For  these  data,  Tj,/Tj  was  of  the  order  of  one. 
The  correlation  based  on  the  measured  shows  the  most  scatter.  The 
correlations  based  on  and  have  about  the  same  scatter,  although  the 
latter  has  the  feature  of  improving  the  correlation  at  low  Mj  (i.  e.  ,  the  data 
of  references  15  and  16). 
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Figure  21  WAKE  TRANSITION  CORRELATION  (BASED  ON  MEASURED 
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Figure  22  WAKE  TRANSITION  CORRELATION  (BASED  ON  COMPUTED 
TRANSITION  DISTANCE) 
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VII.  PREDICTION  OF  TRANSITION 


The  Reg,  Mj  and  Re^  ,  Mj  histories  of  the  wake  of  a  5-foot,  12-degree  cone 
have  been  computed  and  are  shown  in  figures  23  and  24.  *  The  predicted  loci 
of  transition  points  is  remarkably  consistent  for  these  two  correlations. 
(Predictions  on  the  basis  of  Re^^^  indicate  a  long  laminar  wake  persisting  to 
altitudes  lower  than  100,  000  feet).  Presumably,  the  transition  Reynolds  number 
for  the  boundary  layer  will  continue  to  rise  with  Mach  number  inferring  that 
transition  on  the  cone  surface  will  occur  in  the  neighborhood  of  100,000  feet, 
at  which  time  the  near -wake  and  any  portion  of  the  far -wake  still  laminar 
would  be  tripped.  ** 

The  effect  of  increased  vehicle  size  also  has  been  investigated.  As  can  be 
seen  from  figures  21  and  22,  the  transition  point  stays  near  the  neck  at 
115,  000  feet  regardless  of  vehicle  size,  although  as  much  as  300  to  400  feet 
of  laminar  rim  may  occur  behind  a  15 -foot  cone  at  200,  000  feet. 

The  possible  effects  of  chemical  nonequilibrium  on  the  prediction  of  transition 
for  this  vehicle  have  not  been  included;  however,  the  velocity  decay  is  only 
slightly  affected  by  the  temperature  variation  (e,  g.  ,  figure  10). 

It  should  be  noted  for  large  X  the  values  of  and  Reg  become  an  inverse 
function  of  X  .  Hence,  at  some  altitude  above  200,  000  feet,  in  the  present 
example,  the  wake  will  remain  completely  laminar. This  result  is  qualitatively 
consistent  with  conclusions  reacted  by  Lees^  from  energy  considerations. 


*  A  constant  velocity  (22,000  ft/sec)  trajectory  between  200,000  and  100,000  feet  is  assumed. 

••  It  is  anticipated  that  surface  ablation  will  have  two  opposing  effects  on  the  boundary-layer  stability.  The  blowing 
itself  tends  to  be  destabilizing,  but  on  the  other  hand  it  will  enhance  the  stabilizing  effects  of  the  favorable  pressure 
gradient  (induced  by  viscous-inviscid  interaction)  by  further  increasing  the  displacement  thickness. 

***This  altitude  bound  on  the  occurrence  of  transition  will,  of  course,  depend  on  the  size,  shape,  and  velocity  of  the 
vehicle  in  question. 


-47- 


Figure  23  WAKE  TRANSITION  FOR  A  12.DEGREE  CONE  AT  22,000  FT/SEC 

BASED  ON  Reg 
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I  24  WAKE  TRANSITION  FOR  A  12'DEGREE  CONE  AT  22,000  FT/SEC 
BASED  ON  Re 
63-3924  * 
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vm.  NECK  THICKNESS  OF  THE  LAMINAR  WAKE 


As  in  the  previous  discussion  (section  IV)  on  the  momentum  thickness  at  the 
neck,  it  also  may  be  argued  that  the  overall  wake  thickness  at  the  neck  also 
must  vary  as  the  fourth  root  of  Reynolds  number  in  the  axisymmetric  case, 
as  opposed  to  the  square  root  in  the  planar  (semi -infinite)  case.  In  particular, 
if  negligible  mass  addition  to  the  viscous  core  in  the  near-wake  region  is  as¬ 
sumed,  then  a  simple  mass  balance  yields: 
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where 


The  parameter 


_S 

d 


(Rej 


neck 


was  found  to  be  essentially  a  function  solely 


of  Mach  number  for  the  range  of  cone -angle  and  Mach-number  combinations 
encountered  in  the  present  tests,  when  computed  in  the  manner  discussed  in 
in  the  previous  section  on  initial  conditions  (i.  e.  ,  employing  the  values  of 
S  which  give  the  proper  drag  at  the  neck  when  used  with  the  assumed  separa¬ 
tion-point  velocity  profile).  The  values  of  this  parameter  based  on  the  experi¬ 
mentally  measured  neck  thickness  are  shown  in  figure  25  to  compare  remark¬ 
ably  well  with  the  computed  curve.  * 


*The  thickness  S  in  this  figure  is  the  wake  diameter  rather  than  radius. 
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FOR  SHARP 


IX.  CONCLUSIONS 


The  results  of  theoretical  and  experimental  investigations  of  the  flow  in  the  wake 
of  slender  hypersonic  vehicles  have  been  presented.  In  particular,  nonsimilar 
solution  of  the  governing  equations  is  contrasted  with  similar  and  Karman- 
Pohlhausen-type  solutions.  The  similar  solution  is  shown  to  be  inadequate  near 
the  neck  of  the  wake  of  a  slender  body,  although  correct  asymptotically  at  large 
distances.  The  Karman-Pohlhausen  technique  appear  s  to  offer  a  good  compromise 
between  simplicity  and  accuracy  for  slender -body  wakes.  Some  numerical 
examples  for  slender  and  blunt  bodies  are  presented  using  the  nonsimilar  solution. 

The  results  of  a  series  of  ballistic  range  experiments  on  slender  cones  in  the 
Avco  RAD  facilities  are  summarized  and  discussed.  The  preceding  analysis  is 
used  to  compute  appropriate  parameters  for  the  correlation  of  wake  tr«msition 
and  neck  diameter  from  these  data.  Moreover,  predictions  of  wake  transition 
for  a  5-foot,  12-degree  half-angle  cone  at  22,000  ft/sec  are  made,  which  indicate 
transition  moves  from  about  250  feet  behind  the  cone  apex  at  an  adtitude  of 
200,000  feet  to  about  20  feet  at  100,000  feet.  Increasing  the  cone  length  to  15 
feet  is  shown  to  add  about  100  feet  to  the  length  of  the  laminar  trail  at  200,  000 
feet,  and  to  add  about  15  feet  at  100,000  feet.  The  existence  of  an  altitude  above 
which  transition  cannot  occur  (for  a  vehicle  of  given  size  aind  velocity)  adso  is 
indicated  by  the  correlation. 
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